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BENCHMARK PROBLEMS AND SOLUTIONS 
Christopher K.W. Tam 

The Scientific Committee, after careful consideration, adopted six categories of benchmark 
problems for the workshop. These problems do not cover all the important computational issues 
relevant to Computational Aeroacoustics (CAA). The deciding factor to limit tire number of cat- 
egories to six was the amount of effort needed to solve these problems. For reference purpose, the 
benchmark problems are provided below. They are followed by the exact or approximate analyti- 
cal solutions. At present, an exact solution for the Category 6 problem is not available. 


BENCHMARK PROBLEMS 


Category 1 

Problems to test the numerical dispersion and dissipation properties of a computation 
scheme (linear waves). 

Use nondimensional variables with the following scales 

Ax = A r = length scale 

Qoo (ambient sound speed) = velocity scale 
Ax 

= time scale 

Oqo 

Poo = density scale 
Pooaio = pressure scale 


1. Solve the initial value problem 


du du 
dt dx 


t — 0 u = 0.5 exp 


-(‘ n 2 ) (I) 


Give numerical solution at t = 100, 200, 300 and 400 over —20 < x < 450. State the size of 
At used. 


2. Solve the spherical wave problem 


du u du 

^ ^ 7T = ® 

dt r dr 
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over the domain 5 < r < 450, with initial condition t = 0, u = 0. The boundary condition at 
r — 5 is: 

r = 5, u = sin ut 

t ^ x 

(a) w = - 

(b) ^ = f 

Give the numerical solution at t = 100, 200, 300 and 400 for each case. (Do not recast the 
equation in a plane wave form.) State the size of At used. 


Category 2 


Problems to test the nonlinear wave propagation properties of a computational scheme. 
Use dimensionless variables with the following scales 


Ax = length scale 

doc (ambient sound speed) = velocity scale 
= time scale 

®oo 

p 0 o = density scale 
Poodle = pressure scale 


In both problems, the one-dimensional Euler equations are to be solved. 


dp dpu 
dt dx 

( du du\ dp 

^ \ dt U dx ) dx 


dp dp du 

ar + u to +7P & =0 


(You may use an equivalent form of the Euler equations.) 


1. Solve the initial value problem 

t = 0 u = 0.5 exp 
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p = 1 + 


7- 1 


7-1 


7 = 1.4 


Use a computational domain —50 < x < 350. Give the spatial distribution of u, p, and p at 
t = 10, 20, 30, 40, 50, 100, 150, 200 and 300. 


2. Solve the one-dimensional shock tube problem using the following initial conditions 

x < — 2 
—2<x<2 
x > 2 

P = ('YP)'y 7 = 1-4 

Use a computational domain -100 < x < 100. Give the spatial distribution of p, p and u at 
t = 40, 50, 60 and 70. 


t = 0 u = 0 



Category 3 

Problems to test the effectiveness of radiation boundary conditions, inflow and outflow 
boundary conditions and the isotropy property of the computation algorithm. 

Use dimensionless variables with the following scales 

Ax — length scale 

a Q 0 (ambient sound speed) = velocity scale 
Ax 

= time scale 

Qoc 

Poo = density scale 
Po oOm = pressure scale 


In both problems, the linearized two-dimensional Euler equations on a uniform mean flow are 
to be solved. 


dU dE dF 

dt + dx + dy 
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where 


, E = 


' M x p + w 
M x u + p 
M x v 

,M x p + u. 


' MyP + V 
MyU 
MyV + P 
■ MyP + V . 


M x and M y are constant mean flow Mach number in the x and y direction, respectively. 
(You may use an equivalent form of the above equations.) 

Use a computational domain —100 < x < 100, —100 < y < 100 embedded in free space. 
1. Let M x — 0.5, M y = 0. Solve the initial value problem, t = 0. 


p — exp —(In 2) 


x 2 + y 2 


/> = -p f— (In 2) (i±^)| + 0.1 exp 2 


u = 0.04yexp —(In 2) 


(x - 67) 2 + y 2 


v = — 0.04(a: — 67) exp —(In 2) 


{x - 67) 2 + y 2 
25 


Give the distributions of p, p, u and u at t — 30, 40, 50, 60, 70, 80, 100, 200 and 600. 


2. Let M x = M y = 0.5cos( — ). Solve the initial value problem, t = 0. 


p ---- exp —(In 2) 


x 2 + y 2 


, = exp f — (In 2) (*1±^)| +0.1 exp [-(In 2) (» - ~ 


= 0.04(y — 67) exp —(In 2) 


X x ~ 67 ) 2 + (y - 67 ) 2 


v = — 0.04(x — 67) exp —(In 2) 


.(*-6 7) 2 +(y-67) : 
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Note: The mean flow is in the direction of the diagonal of the computational domain. 
Give the distributions of p, p, u and u at t — 60, 70, 80, 90, 100, 200, 600 and 1000. 




Category 4 

Problems to test the effectiveness of wall boundary conditions. 

Use dimensionless variables with the following scales 

Ax = Ar = length scale 

Ooo (ambient sound speed) = velocity scale 
Ax 

= time scale 

a oo 

Poo = density scale 
Poodle ~ pressure scale 


1 . 


Reflection of an acoustic pulse off a wall in the presence of a uniform flow in semi-infinite 
space. 
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Use a computational domain —100 < x < 100, 0 < y < 200. The wall is at y = 0. The 
linearized Euler equation in two dimensions are 





' Mp + u ' 


~v‘ 

d 

u 

d 

^ dx 

Mu + p 

d 

+ dy 

0 

dt 

V 

-P- 

Mv 

. Mp -f u . 

P 

.v . 


where M = 0.5. The initial condition is 


t — 0, u = v = 0 

x 2 4- (y - 25) 2 ' 


p — p = exp < —(In 2) 
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} 


Give the pressure field at t = 15, 30, 45, 60, 75, 100 and 150. 


2. Acoustic radiation from an oscillating circular piston in a wall 



7rt 

Radius of piston = 10. Velocity of piston it = 10~ 4 sin( — ). Use a computational domain 

5 

0 < x < 100, 0 < r < 100. The wall and the piston are at x = 0. The cylindrical coordi- 
nate system is centered at the center of the piston. With axisymmetry, the linearized Euler 


i t; 
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equations are 


d_ 

dt 


' p' 


‘ v' 


r v 
r 


~u~ 

u 

d 

-f- — — 

0 

+ 

0 

d 

+ — — 

p 

V 

dr 

P 

0 

dx 

0 

-P- 


. v . 


V 

L r J 


,u_ 


The initial conditions are: 


t = 0 


p = u — v=p = 0 


= 0 


Give the time harmonic pressure distribution at the beginning, 
piston oscillation. 


\ and 7 ofa period of 
2 4 


Category 5 


Problem to test the suitability of a numerical scheme for direct numerical simulation of very 
small amplitude acoustic waves superimposed on a non-uniform mean flows in a semi-infinite 
duct. 

Use nondimensional variables with the following scales 


Ax = length scale 

a 0 o (sound speed far upstream) = velocity scale 

j/\ T 

= time scale 

Ooo 

Poo (density of gas upstream) = density scale 
pooO^oo = pressure scale 


A small amplitude sound wave is incident on a convergent -divergent nozzle as shown 


M=0.5 _ . Sound Wave 

— ~ / VU— 



Supersonic 



Figure for Category 5 Problem 
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Use a computational domain —200 < x < 80. The area of the nozzle is given by 


A(x) = < 


f 134 

117-17cos(^) -100<x<19 

(97.2 + 0.3x 


x < -100 
-100 < x 
19 < x < 80 


The quasi-one-dimensional unsteady flow equations are 

dp A dpuA 


+ 


= 0 


dt dx 
( du du\ dp 

p \dt JrU Tx) + Tx =a 

Far upstream x < —200, there is an incoming acoustic wave. Together with the steady inflow, the 
velocity, pressure and density are given by 


u 


'M‘ 


V 



/ \ n 

p 


1 

+ 

1 

e sin 

U) ( 

X i\ 
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yl + M )_ 

_p_ 


1 


_1 





Take 7 = 1.4, M = 0.5, e = 10 6 , u = 0.17T, calculate the transmitted sound wave at the nozzle 
exit. Give p(t) — p over a time period; p is the time averaged pressure. 


Category 6 


Problems to test the ability of a numerical scheme to calculate aeroacoustic source. 


Use dimensionless variables with the following scales 


Ax = length scale 

aoo (ambient sound speed) = velocity scale 

Ax . , 

= time scale 

Qoo 

p a o = density scale 
PooO^x, = pressure scale 


1. Sound generation by gust-blade interaction (two-dimensional) 


8 
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Figure for Category 6 Problem 

Use a computational domain —100 < x < 100, —100 < y < 100. The blade is a flat plate of 
length L (L = 30) lying along the x-axis centered at the origin. There is a Mach 0.5 uniform 
mean flow in the x-direction. The mean flow carries a gust with velocity component in the 
y-direction given by 

” = 01sin [f Gt;-')]’ M “ = 0 - 5 

The linearized Euler equations are 

'pi rMoop + ui vv 

dud Moou +p d 0 
dt v dx MqqV Qy p 

■ p j L MooP + u J iv. 

Determine the intensity of radiated sound, p 2 , along the coordinate lines x = ±95 and y = 
±95. 

SOLUTIONS 
Category 1 

Problem 1. The solution is 

u(x, t) — 0.5 exp —(in 2) 

Problem 2. The solution is 




u ( r >0 | |[sina;(< - r ± 5)], 


r > t + 5 
r < t + 5. 
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Category 2 


Problem 1. An approximate solution can be found by using the simple wave assumptions (Chap- 
ter 6, G.B. Whitham, “Linear and nonlinear waves”). These assumptions are 

1. The flow is isentropic 

2. The Riemann invariant — u = — j-, which starts from the uniform region ahead of the 
pulse, is valid everywhere. 

With these assumptions, the Euler equations reduces to the nonlinear simple wave equation 

du ( 7 + 1 \ du 

dt + { 1 + ~ u ) di - 0 

This quasi-linear first-order equation can be solved by the method of characteristics. For the 
given initial conditions, a shock will form at the front of the pulse as the disturbance propagates 
to the right. The location of the shock may be found approximately by the use of Whitham’s 
equal area rule. 

Problem 2. The standard shock tube solution is a good approximate solution. The standard so- 
lution is available in most books on gas dynamics. 


Category 3 


Problem 1. Let aq = a 2 = M = 0.5, r] = [(x — Mt ) 2 + y 2 ]? 
The solution is 


(x - Mt) 


sm((t)Ji((r])£ d£ + 0. 

i *7 J 


-Q2[(i- 67 — Mt) s +5i 2 


2a iy 


OO 

J e^T gin (£<) (£>?)£ d£ - 0.04(x - 67 - Mt)e- a ^ x - 67 ~ Mt ^ +y ^ 


OO 2 

P = g— J cos(£f)J 0 (£»?)£ <*£ 


/ ) = p + 0.1e-“^- 67 - Af<)2+ ^ 

where Jq( ) and Ji( ) are Bessel functions of order 0 and 1. 
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(Reference: C.K.W . Tam and J.C. Webb, “Dispersion-Relation-Preserving finite difference 
schemes for computational acoustics,” J. Computational Phys., Vol. 107, pp. 262-281, 1993.) 

Problem 2. The solution can be obtained from that of problem 1 by a coordinate transformation. 


Category 4 

Problem 1. Let a = rj = [(x - Mt ) 2 + (y - 25) 2 ] i, £ = [(*- Mt) 2 + (y + 25) 2 ] 
The solution is 


u — 


( x — Mt) [_&!., 


2ar) 


J e sin(^)Ji(^X + — 2 ~^ J e ^ sin(#)-M£C)£ # 


y = 


(y - 25 ) 

2arj 


00 00 

J e-^sin + J 


e ^sin 


=p= h] 


. lL 


e *« cos(^)[J 0 (^) + JoUC)]C 


Problem 2. Let e = 10 4 , R — 10, u: = f. 
The solution is 


p = p = Re 


eRlo 


OO 

I 




(£ 2 — ^> 2 ) 2 


u = Im 


u = Im 


-eR 


OO 

J M(R)M(r)> 


— (£ 2 — w 2 ) 2 x — iuit 


-eR 


00 

J 


Jim 




r £Ji(£r)e 




— (£ 2 — u/ 2 ) 2 i — iu>< 




where Re[ ] = the real part of and Im[ ] = the imaginary part of. 
Note: (£ 2 -w 2 )i = -i'|£ 2 -u; 2 |5 for £ < u. 
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Category 5 


A fairly accurate solution of this problem can be found by first determining the governing 
equations for the amplitude functions of the time-periodic disturbances inside the nozzle. These 
equations are ordinary differential equations but with variable coefficients. They can be inte- 
grated numerically. 

Let the solution be separated into a mean and a time-periodic part in the form 


' p 


p 

f 

'p(x)' 

e-'“‘ 

u 


u 

+ Re < 

u{x) 

,p_ 


.p. 

1 

_p( x )_ 

J 


The physical quantities of the mean flow at the nozzle throat will be denoted by a subscript *. 
With the area ratio A*/A\ known, where A\ is the area of the uniform duct, p* is first found by 
solving the equation 


7 +i / A±\ + _£_ 

P * 1 A x ) 7-1 7-1 


— rt' 1 


The other variables at the nozzle throat are given by 

1 
7 


p* = -pi, 


xni 

U* = P* 2 


The mean flow solution is 


puA = p*u*A* 


P ( P x 7 


P * 


P* 


_ 7 __ P*_ 797—1 _ 


= — + 


7 P* 


2 7 — 1 pl P 2 7 - 1 p* 

The linearized governing equations for the amplitude functions p, u and p are 

1 


du 

dx 

dp 

dx 


(p v? -7 p) 

1 

(pu 2 -7 p) 


du ^ (. __du dp 7 pdA\^ (. d(u A) 


A dx 


du ^ ( ; • du dp jppudA 

1PU li P+ t-"TP?+7PA ^ - P* Tx - — Tx 


U 


/. 7 pu d(u A) 

+ s~IP 


( 1 ) 


dp p du (. d(uA)\ 1 ^ 1 d(p A) ^ 

—— = — — — j - I iuiA — — ) zt~t p — r , n 

dx u dx \ dx J u A u A dx 
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In the uniform region of the duct, the solution of (1) that matches the incoming acoustic 
wave is 


r 

p 


V 


1 ‘ 

u 

= — is 

1 

tU5X 

-j. c 

-1 

.p. 




1 


( 2 ) 


In (2), the second term represents the reflected acoustic wave. The unknown amplitude c is to be 
determined later. 

Equations (1) have a regular singular point at the nozzle throat (x = 0). Near the throat, 
there are two non-singular series solutions. The first two terms of these solutions are 


P = 


xujp * 


u 


if 1 

*47+1 )A. 


(—) 

V dx 2 )* 


2/5* 


U, 


Uo + 


-UjJ 


u 


ar i 

* 47 +l)A,t dx 2 


f£4). 15 


+ 


7 


ut 


Po + PlX + 


U = Uo + UiX + . 


( 3 ) 


P = Po + P\X + ... 

where u 0 and po are arbitrary constants, px, ux and px are functions of u 0 and p 0 . 

A numerical solution of (1) can be constructed by starting the solution slightly upstream of 
the nozzle throat at x = —8 ( 8 « 1) using (3) as the starting solution. (For small 5, the terms 
of the series involving 6 and powers of 8 may be neglected.) The numerical integration proceeds 
upstream until the uniform duct region is reached. At this point, the numerical solution must 
match solution (2). This provides three algebraic equations for the three unknowns po, and c. 
Once these constants are found, the solution upstream of the nozzle throat is known. 

For the solution downstream of the nozzle throat one can start integrating (1) numerically at 
a point just downstream, say at x = 8. Again (3) is used as the starting solution. The numerical 
integration proceeds downstream until the nozzle exit is reached. With p 0 , uo already found, the 
amplitude functions are now completely determined along the entire length of the nozzle. 


Christopher K.W. Tam 
Department of Mathematics 
Florida State University 
Tallahassee, FL 32306-3027 
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